The Symmetry Group of Gaussian States in 

l2(M") 



K. R. Parthasarathy 



Abstract This is a continuation of the expository article f4\ with some new re- 
marks. Let S„ denote the set of all Gaussian states in the complex Hilbert space 
L?{W), Kf, the convex set of all momentum and position covariance matrices of or- 
der 2« in Gaussian states and let be the group of all unitary operators in L"(R") 
conjugations by which leave S„ invariant. Here we prove the following results. K„ is 
a closed convex set for which a matrix S is an extreme point if and only if S = jL^L 
for some L in the symplectic group Sp{2n,M.). Every element in K„ is of the form 
^{L^ L + M'^ M) for some L,M in Sp{2n,R). Every Gaussian state in L^{R") can 
be purified to a Gaussian state in L^(R^"). Any element U in the group is of 
the form U = XW {(X)r (L) where A is a complex scalar of modulus unity, CC G C", 
L G 5'p(2«,R), W{(X) is the Weyl operator corresponding to (X and F(L) is a uni- 
tary operator which implements the Bogolioubov automorphism of the Lie alge- 
bra generated by the canonical momentum and position observables induced by the 
symplectic linear transformation L. 
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1 Introduction 



In in we defined a quantum Gaussian state in L^(IR") as a state in which ev- 
ery real linear combination of the canonical momentum and position observables 
p\,P2, ■ ■ ■ iPiu <Jii<j2, ■■■,qn has a normal distribution on the real line. Such a state 
is uniquely determined by the expectation values of pi,p2, ■■■ ,Pn', qi,q2, ■■■ ,qn and 
their covariance matrix of order 2«. A real positive definite matrix S of order 2n is 
the covariance matrix of the observables p\,p2, ■ ■ ■ ,Pn', qi,q2, ■ ■ ■ ,qn if and only if 
the matrix inequality 

2S-iJ>0 (1.1) 



holds where 



-/ 
/ 



(1.2) 



the right hand side being expressed in block notation with and / being of order 
nxn. We denote by K„ the set of all possible covariance matrices of the momentum 
and position observables in Gaussian states so that 

A'„ = {5'|5' is a real symmetric matrix of order 2« and 25 — //> 0} . (1.3) 

Clearly, K„ is a closed convex set. Here we shall show that S is an extreme point of 
K„ if and only if 5 = ^L^L for some matrix L in the real symplectic matrix group 

Sp{2n,R) = {L\L^JL^J} (1.4) 

with the superfix T indicating transpose. Furthermore, it turns out that every element 
S in K„ can be expressed as 

s=^{l^l+m'^m) 

for some L,M € Sp{2n,'R). This, in turn implies that any Gaussian state in L^(M") 
can be purified to a pure Gaussian state in L^(M^"). 

Let a e (ai , a2, . . . , a„y e C", L = ((£,;■)) G Sp{2n,R) and let aj ^xj + iyj with 
Xj,yj £ M. Define a new set of momentum and position observables p\,p'2, ■ ■ ■ ,p'n', 

q'i,q'2i---^i'i^y 

n 

p'i = L {^ijiPj-'''j) + ^in+Mj~yj)] , 

n 

^' = L {^n+ijiPj-^i)+^"+i"+jili~yi)] ^ 

for 1 < i < n. Here one takes linear combinations and their respective closures to 
obtain p-,^- as selfadjoint operator observables. Then p[, p^,- ■ ■ , p',{, ^'■,^'2, . . . 
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obey the canonical commutation relations and thanks to the Stone-von Neumann 
uniqueness theorem there exists a unitary operator F( a, L) satisfying 

p'i = r{a,L)piria,L)\ 
q'i^r{a,L) qiria,L)^ 

for all 1 < / < n. Furthermore, such a F(a,L) is unique upto a scalar multiple of 
modulus unity. The correspondence (a,L) — > r{a,L) is a projective unitary and 
irreducible representation of the semidirect product group C" (S) Sp{2n,'R). Here 
any element L of Sp{2n,'R) acts on C" real-linearly preserving the imaginary part 
of the scalar product. The operator r{a,L) can be expressed as the product of 
W{a) = r(a,/) and r(L) = r(0,L). Conjugations by W{a) implement transla- 
tions of Pj,qj by scalars whereas conjugations by F(L) implement symplectic lin- 
ear transformations by elements of 5'/9(2«,R), which are the so-called Bogolioubov 
automorphisms of canonical commutation relations. In the last section we show that 
every unitary operator U in (M" ) , with the property that UpU^ is a Gaussian state 
whenever p is a Gaussian state, has the form U = XW {a)r (L) for some scalar A of 
modulus unity, a vector a in C" and a matrix L in the group Sp{2n,R). 

The following two natural problems that arise in the context of our note seem to 
be open. What is the most general unitary operator U in L?{W) with the property 
that whenever is a pure Gaussian state so is Secondly, what is the most 

general trace-preserving and completely positive linear map A on the ideal of trace- 
class operators on L^(M") with the property that A (p) is a Gaussian state whenever 
p is a Gaussian state? 



2 Exponential vectors, Weyl operators, second quantization and 
the quantum Fourier transform 

For any z = (zi,Z2, ■ ■ ■ ,Zn)^ in C" define the associated exponential vector e{z) in 
by 

e{z){x) = (2;r)-"/4exp £ [zjXj \z) \x)). (2.1) 

Writing scalar products in the Dirac notation we have 

(.(z)k(z')) =exp(z|z') 

n 

= exp^Z;Z^-. (2.2) 

7=1 

The exponential vectors constitute a linearly independent and total set in the Hilbert 
space l? (R" ) . If t/ is a unitary matrix of order n then there exists a unique unitary 
T{U) inL2(R«) satisfying 
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r{U)\e{z)) ^\eiUz)) VzeC". (2.3) 

The operator r{U) is called the second quantization of U. For any two unitary 
matrices U,V in the unitary group ^ (n) one has 

r{u)r{v)=r{uv). 

The correspondence U r(U) is a strongly continuous unitary representation of 
the group («) of all unitary matrices of order n. 

For any a £ C" there is a unique unitary operator W(oc) in L'(K") satisfying 

lV(a) |e(z)) |e(z + a)) V zeC". (2.4) 

For any Of , j3 in C" one has 

W(a)W(/3) = W(a + /3). (2.5) 

The correspondence a W{a) is a projective unitary and irreducible represen- 
tation of the additive group C". The operator W{(X) is called the Weyl operator 
associated with a. As a consequence of ( 12.5b it follows that the map t W{t(X), 
f € R is a strongly continuous one parameter unitary group admitting a selfadjoint 
Stone generator p{(X) such that 

w{ta) = e-"p^"^ V aeC". (2.6) 

Writing = (0, 0, . . . , 0, 1 , 0, . . . , 0)^ with 1 in the ;-th position, 

pj=2-^p{ej), qj = -2-^p{iej) (2.7) 

(2.8) 

one obtains a realization of the momentum and position observables pj,qj,l < i <n 
obeying the canonical commutation relations (CCR) 

[Pi^Pj] = 0, [qi,qj] = 0, [q,-,p.s] = iSrs 
and the adjoint pairs aj, a^j of annihilation and creation operators satisfying 

[a,,fl;]=0, [«;,«]] = 5,7 
in appropriate domains. If we write 

p'j = 2-^Pj, q'j^l^qj 

one obtains the canonical Schrodinger pairs of momentum and position observables 
in the form 
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{p'jV) W = j S'/^^' ^^^"^^ ^''^ ^^^^ 

in appropriate domains. We refer to iQ for more details. 

We now introduce the sympletic group Sp{2n,R.) of real matrices of order 2n 
satisfying ilAi . Any element of this group is called a symplectic matrix. As de- 
scribed in IT], im, for any symplectic matrix L there exists a unitary operator F(L) 
satisfying 

r{L)w{a)r{Ly ^w{La) V aeC" (2.9) 



where 



Re La 
Im La 



= L 



Re a 
Im a 



(2.10) 



Whenever the symplectic matrix L is also a real orthogonal matrix then Lisa unitary 
matrix and F(L) coincides with the second quantization F(L) of L. Conversely, if 
L'^ is a unitary matrix of order n, Lu is the matrix satisfying 



Lu 



Re U{x + iy) 
Im L'^ (x + iy) 



then Lu is a symplectic and real orthogonal matrix of order 2n and r{Lu) = r{U). 
Equations ( 12.91 ) and ( 12.6b imply that F(L) implements the Bogolioubov automor- 
phism determined by the symplectic matrix L through conjugation. 

For any state p in L^(M") its quantum Fourier transform p is defined to be the 
complex-valued function on C" given by 



p{a) = TrpW{a), aeC". (2.11) 

In ||4] we have described a necessary and sufficient condition for a complex-valued 
function / on C" to be the quantum Fourier transform of a state in L" (M" ) . Here 
we shall briefly describe an inversion formula for reconstructing p from p . To this 
end we first observe that (12.11b is well defined whenever p is any trace-class op- 
erator in L^(R). Denote by and ^2 respectively the ideals of trace-class and 
Hilbert- Schmidt operators in L^(R"). Then JFi C ^2 and ^2 is a Hilbert space 
with the inner product {A\B) = TrA^B. There is a natural isomorphism between ^2 
and L^(M") (X)L^(M"), which can, in turn, be identified with the Hilbert space of 
square integrable functions of two variables x,y in R". We denote this isomorphism 
by so that J^{A){x,y) is a square integrable function of {x,y) for any A € ^2 and 

^(|e(«))(.(v)|)(x,3') -.(«)W.(v)W (2.12) 
for all u. V C C". v denoting (v"! , v'2, . . . , v„). From ( 12.4b and (12. lib we have 



\e{u)){e{v)\r{a) = {e{v)\W{a)\e{u)) 



exp 



(ecu) 
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Substituting a ~x + iy and using (12.1b . the equation above, after some algebra, can 
be expressed as 

i\eiu){eiv)\r{x + iy) = (2;r)"/2e(«')(V2;c)e(v')(V2y) (2.13) 

where 





= U 


u 


v' 




V 



U = 2-1/2 



/ / 
// il 



(2.14) 



LetDe, > denote the unitary dilation operator in L'^{R") igi L^{R") defined by 



iDef){x,y)^e"f{ex,9y). 
Then (12.13b can be expressed as 



(2.15) 



{\e{u)){eiv\nx + iy) = 7t"/^{D^,riU)e{u®v)}{x,y) 

where F(C/) is the second quantization operator in L?(S?") associated with the uni- 
tary matrix U in (12.14b of order 2«. Since exponential vectors are total and and 
r{U) are unitary we can express the quantum Fourier transform p ^ p{x + iy) as 

p = n"/^D^r{U)J^{p). (2.16) 
In particular, p{x + iy) is a square integrable function of {x,y) G M." x W and 

p = 7t-"/^ J^-^ r{U^)D2-,/2P (2.17) 

is the required inversion formula for the quantum Fourier transform. It is a curious 
but an elementary fact that the eigenvalues of U in ( 12. 14b are all 12* roots of unity 
and hence the unitary operators r{U) and r{U^) appearing in (12.16b and ( 12.17b 
have their 12-th powers equal to identity. This may be viewed as a quantum analogue 
of the classical fact that the 4-th power of the unitary Fourier transform in L^(M.") is 
equal to identity. 



3 Gaussian states and their covariance matrices 



We begin by choosing and fixing the canonical momentum and position observables 
pi,/92, ••■ ?i,?2, •••,?)! as in equation (12.7b in terms of the Weyl operators. They 

obey the CCR. The closure of any real linear combination of the form ^ {xjpj — 

7=1 

yjqj) is selfadjoint and we denote the resulting observable by the same symbol. As 
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in El, for a = («! , a2, • • • , 0C„)^, aj — xj + iyj with xj, yj <E M, the Weyl operator 
W{a) defined in Section|2]can be expressed as 

-iV2 I, {xjPj-yjCij) 
W{a)=e J=' . (3.1) 

Sometimes it is useful to express W{0l) in terms of the annihilation and creation 
operators defined by (12.8b : 

W{a)^e'=' ■ (3.2) 

where the linear combination in the exponent is the closed version. A state p in 

n 

L (R) is said to be Gaussian if every observable of the form ^ {xjPj —yjqj) has a 

7=1 

normal distribution on the real line in the state p for Xj,yj e M. From ||4| we have 
the following theorem. 

Theorem 1. A state p in L^(M") is Gaussian if and only if its quantum Fourier 
transform p is given by 



p{a) = TrpW(a) 




(3.3) 



for every (X = x + iy, x,y Cz M" where i, tn are vectors in W and S is a real positive 
definite matrix of order 2n satisfying the matrix inequality 2S — iJ > 0, with J as in 
(fOI) . 

Proof. We refer to the proof of Theorem 3 . 1 in Q . □ 

We remark that £, m and S in ( 13.31 ) are defined by the equations 

n 

e^x - tri^y = Tr p £ {xjpj - yjqj) 
j=i 

{x^y)S (^^) = Tr p X2 - (Tr pX)2,X = ti^jPj - yjqj). 

It is clear that ij is the expectation value of pj, mj is the expectation value of qj 
and S is the covariance matrix of pi,p2, ■ ■ ■ ,pn', ~qi,—q2,---,—qn in the state p 
defined by ( I3.3l l. By a slight abuse of language we call S the covariance matrix of 
the Gaussian state p . All such Gaussian covariance matrices constitute the convex 
set Kn defined akeady in (11.3b . We shall now investigate some properties of this 
convex set. 

Proposition 1 (Williamson's normal form |[1]). Let A be any real strictly positive 
definite matrix of order 2n. Then there exists a unique diagonal matrix D of order n 
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with diagonal entries d\>d2>->dn>0 anda symplectic matrix M in Sp(2n,'R.) 
such that 



M' 



D 
D 



M. 



(3.4) 



Proof. Define 



where J is given by ( 11.2b . Then B is a real skew symmetric matrix of full rank. Hence 
its eigenvalues, inclusive of multiplicity, can be arranged as ±idi,±id2, ■ ■ ■ ,iid„ 
where di > d2 > ■ ■ ■ > d„ > 0. Define D = diag((ii ,c/2, • • . ,dn), i.e., the diagonal 
matrix with dj as the ii-th entry for 1 < i < n. Then there exists a real orthogonal 
matrix F of order 2n such that 



Define 



L^A'i'-r 



Then LfjL^J and 



LAL' 



-D 
D 





D-i/2 



D 
D 



Putting M = (L" ' )^ we obtain (Il4l l. 

To prove the uniqueness of D, suppose D' = diag{d[,d'2, ■ ■ ■ ,d'^) with d\ > c/t — 
■ ■ ■ ^d'„ > and M' is another symplectic matrix of order 2n such that 



A=M' 



D Q' 




D' 


" 


m = m'^ 


D 





D' 



M'. 



Putting = MM we get a symplectic such that 



D 
D 



Substituting A^' = JN'^J' ' we get 



A^" 



D 
D 



A^ = 



A': 



D' 
D' 



D' 

-D' 



Identifying the eigenvalues on both sides we get D ^ D' □ 

Theorem 2. A real positive definite matrix S is in K„ if and only if there exists a di- 
agonal matrix D = diag{di ,d2, ■ ■ ■ jdn) with di > d2> ■ ■ ■ > d„ > j and a symplectic 
matrix M G Sp{2n, R) such that 
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S = M' 



D 
D 



M. 



In particular. 



AetS = Y\d]>A-'\ 



(3.5) 



(3.6) 



Proof. Let 5 be a real strictly positive definite matrix in Kn . From ( ll.3l l we have 
S> and therefore, for any L G 5'/?(2«,R), 



U SL > -J. 



(3.7) 



D 
D 



D 1/ 



Using Proposition[T]choose L so that 

where D = diag(c/i,c/2, • • ■,dn), di > d2> ■ ■ ■ > dn > 0. Now ( 13.7b implies 

>0. 

The minor of second order in the left hand side arising from the jj, jn + j, n + yy, 
n + jn + j entries is of? — | > 0. Choosing L = ' we obtain (13. 5t and (13. 6t . Now 
we drop the assumption of strict positive definiteness on S. From the definition of 
K„ in ( 11.3b it follows that for any S E K„ one has S + el € K„ for every e > 0. Since 
5 + £/ is strictly positive definite det5 + el > 4"" V £ > 0. Letting e we see 
that ( 13.61 ) holds and S is strictly positive definite. 

To prove the converse, consider an arbitrary diagonal matrix D = diag(c/i,c/2, • • • ,dn) 
with di > d2> ■ ■ ■ >d„ > J. Clearly 



D 0' 




-/" 


D 


— i 


/ 



>o, 



and hence for any M £ Sp{2n, 
In other words. 



D 
D 



M-i 



-/ 
/ 



>0. 



M' 



D 
D 



MeK„ MeSp{2n,\ 



Finally, the uniqueness of the parameters di > d2> ■ ■ ■ > d„ > j in the theorem is 
a consequence of Proposition[T] □ 
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We now prove an elementary lemma on diagonal matrices before the statement 
of our next result on the convex set K„. 

Lemma 1. Let D > I be a positive diagonal matrix of order n. Then there exist 
positive diagonal matrices Di,D2 such that 

D=^-{D,+D2) = \{D-' +D^'). 

Proof. We write D2 ~ D\X and solve for D\ and X so that 
2D^Dy{I + X)=D^\l + X-^), 
Dy and X being diagonal. Eliminating D\ we get the equation 

which reduces to the quadratic equation 

X^ + {2-AD^)X + I = Q. 
Solving for X we do get a positive diagonal matrix solution 

X=/ + 2(d2- 1)+2D(d2 -7)1/2. 

Writing 

Di=2D{I + xy\ D2=DiX 
we get Di,D2 satisfying the required property. □ 

Theorem 3. A real positive definite matrix S of order 2n belongs to K„ if and only if 
there exist symplectic matrices L,M such that 

S=-(L^L + M^M). 
4 

Furthermore, S is an extreme point ofK„ if and only ifS = jL^Lfor some symplectic 
matrix L. 

Proof. Let S K,,. By Theorem|2]we express S as 



S=N^ 



D 
D 



N (3.8) 



where is symplectic and D = diag(iii,ii2, ... ,d„), di > d2> ■ ■ ■ > d„ > j. Thus 
2D > I and by Lemma[T]there exist diagonal matrices Di > 0, D2 > such that 

2D=i(Di+D2) = i(Dr'+£>2'). 

We rewrite (13.8b as 
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tN' 





Di 




■£>2 




( 


Dj-^ 


+ 





) 



\N. 



Putting 



we have 



D 



1/2 
1 



D 



-1/2 



1/2 
2 

d; 



-1/2 



Since 



D 



1/2 



D; 



-1/2 



S ^ -{UL + M^M). 



,i= 1,2 are symplectic it follows that L and M are symplectic. 



7 0' 




-/" 


/ 




/ 



This proves the only if part of the first half of the theorem. 
Since 

" f nl fn t' 

>0 

multiplication by on the left and L on the right shows that L^L — iJ >Q for any 
symplectic L. Hence jL^L G K„ L & 5p(2«,R). Since K„ is convex, |(L^L + 
M^M) E K„ , completing the proof of the first part. 

The first part also shows that for an element S of K„ to be extremal it is necessary 
that S = \L^L for some symplectic L. To prove sufficiency, suppose there exist 
L e 5/9(2«,M), 5i,52 € Kn such that 

iL^L=i(5i+52). 
By the first part of the theorem there exist Lj e 5'p(2«, M) such that 



(3.9) 



where Si = \{L\Lx +^2^2), ^2 = ^Ci-si-a +^4^4)- Left multiplication by (L^) 
and right multiplication by L ' on both sides of ( 13.91 ) yields 

1 4 



r^-I 



(3.10) 



7=1 



where 



Mi 



Each My is symplectic and positive definite. Multiplying by J on both sides of (I3.1QI I 
we get 
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Thus 



which impHes 




or 



Mj=I V 1 < y < 4 



Thus 



L^Lj^l'l V j 



and =52. This completes the proof of sufficiency. □ 

Corollary 1. Let Si,S2 be extreme points of Kn satisfying the inequality Sy > Si- 
Then Si = Si- 

Proof. By Theorem [3] there exist L,- G 5'p(2«,M) such that 5; = ^LJ Li, i = 1,2. 
Note that M = LiL^^ is symplectic and the fact that > 52 can be expressed as 
M^M < L Thus the eigenvalues of M^M lie in the interval (0, 1] but their product 
is equal to (detM)^ = 1. This is possible only if all the eigenvalues are unity, i.e., 
M^M = L This at once impHes L\Li ~ L2L2- □ 

Using the Williamson's normal form of the covariance matrix and the transforma- 
tion properties of Gaussian states in Section|3]of we shall now derive a formula 
for the density operator of a general Gaussian state. As in |4j denote by Pg{£,m,S) 
the Gaussian state in L^ (R" ) with the quantum Fourier transform 



where £, m E W and S has the Williamson's normal form 



with M G Sp{2n,R), D = dmg{di,d2, - - - ,dn), di > d2> ■ ■ ■ > d„> j - From Corol- 
lary 3.3 of Q we have 
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V2 



t 



PgH,m,S)W 



m + ii 

~7T 



and Corollary 3.5 of ID implies 

p,(o,o,5) = r(M)-V«(o,o, 



D 
D 



:P,(0,0,5) 



r(M). 



Since 



D 
D 



is a diagonal covariance matrix 



P. 0,0, 



D 
D 



7=1 



where the y'-th component in the right hand side is the Gaussian state in L (R) with 
means and covariance matrix djl2,l2 denoting the identity matrix of order 2. If 



dj = 5 we have 



p,(0,0,i/2) = |e(0))(e(0)| inL2 



If dj > 1/2, writing dj ~ \coihjSj, one has 

p,iOAdjl2)^{l-e-'')e-^^"'" 

with Qjo' ,p,q denoting the operator of annihilation, creation, momentum and posi- 
tion respectively in L^(]R). We now identify L^{M.") and L^{M.f'" and combine the 
reductions done above to conclude the following: 

Theorem 4. Let Pg{i,m,S) be the Gaussian state in L^(M") with mean momentum 
and position vectors m respectively and covariance matrix S with Williamson 's 
normal form 



D 
D 



M, Me5/?(2«,R), 



D = diag(£/i,c/2, ■ • -^d,,), d\ > d2> ■ ■ ■ > dm> d^+i = <im+2 = • • • = "in = 5, dj 
jCOth\sj, 1 < j <m, Sj > 0. Then 



m + // 

p,{e,m,S) = W{—j^)r{M)-' Y\{l-e-^') X 
V2 j=i 

e-^U'i"]"i(^{\e{Q)) X {e{0)\f'"""r{M)W{ 



m + ii 



(3.11) 



where H^(-) denotes Weyl operator, r(M) is the unitary operator implementing the 
Bogolioubov automorphism of CCR corresponding to the symplectic linear trans- 
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formation M and |e(0)) denotes the exponential vector corresponding to in any 
copyofL^iR). 

Proof. Immediate from the discussion preceding the statement of the theorem. □ 

Corollary 2. The wave function of the most general pure Gaussian state inL^{W) 
is of the form 

|v/)=w(a)r(c/)|a,)kAj---k;i„> 



where 



ei{x) = (27r)"'/'^A"'''2exp-4"'A"V, jce R,A > 0, 



U e C", U is a unitary matrix of order n, W{Cl) is the Weyl operator associated 
with 0£, r{U) is the second quantization unitary operator associated with U and 
•^jj 1 ^ 7 ^ n are positive scalars. 

Proof. Since the number operator a^a has spectrum {0, 1,2, . . .} it follows from 
Theorem |4] that Pg{£,m,S) is pure if and only if m = in ( 13.1 It . This impUes that 
the corresponding wave function | can be expressed as 



\xif)=W{a)r{M)-'{\e{0))) 



(3.12) 



where M e Sp{2n,M.) and a = The covariance matrix of this pure Gaussian 
state is ^M^M. The symplectic matrix M has the decomposition yj 



M = Vi 



D 




V2 



where Vi and V2 are real orthogonal as well as symplectic and D is a positive diag- 
onal matrix of order « . Thus 









V2 



where 



D 




V2. 



Since the covariance matrix of | \f/) in ( 13.12b can also be written as jN'^N, modulo 
a scalar multiple of modulus unity | xj/) can also be expressed as 



\xtf)=w{a)r{V2r'r 



D-i 
D 



(3.13) 



If U is the complex unitary matrix of order n satisfying 
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U{x + iy) =x' + iy', 
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x' 


= vl 


X 


y. 







V x,y e 



and D ' = diag(Ai,A2, . . . ,A„) we can express (13.13b as 



\yf)=w{a)r{u) 



Xj 
Ari 



em 



= W{a)r{U)\e^^)\e^^)---\e^J 



where we have identified L^{R") with L^{R)'^" . 

We conclude this section with a result on the purification of Gaussian states. □ 

Theorems. Let p be a mixed Gaussian state in L^{W). Then there exists a pure 
Gaussian state 1 1/) in {W) ® {W) such that 

p =Tr2 f/ \W){WW'' 

for some unitary operator U in L? (R" )®L?{W^) with Tr2 denoting the relative trace 
over the second copy ofL?{W). 

Proof. First we remark that by a Gaussian state in L^(M") (X)L^(M") we mean it 
by the canonical identification of this product Hilbert space with L^(M^"). Let p = 
Pg{£,m,S) where by Theorem[3]we can express 

5=^(l[Li+4L2), LuL2£Sp{2n,R). 
Now consider the pure Gaussian states, 

\WL.)=r{Lt)-' \eiO)), i=l,2 
in L?{W) and the second quantization unitary operator i7) satisfying 

+ V M 



Fo e{u®v) = e 



Vu,veC" 



V2 V2 

in l2(R2") identified with L^iE.") (E}L^{M."), so that 

e{u®v) =e{u)(E)e{v). 
Then by Proposition 3.11 of |4j we have 

Tr2 To (I v/L, ) (V^L J ® I {Wl, |) r^,^ = P, (0, 0, S) . 
If a = ^ we have 

\/2 
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Wia)pg(0,0,S)W{a)^ = Pgie,m,S). 

Putting 

u = (w(a) ®/)ro {r{Li)-'(^r{L2)-') 

we get 

Pg{e,m,S) = Tt2U \e{0)(g,e{0)){e{0)^e{0)\ 
where |e(0)) is the exponential vector in L^(R"). □ 



4 The symmetry group of the set of Gaussian states 

Let Sn denote the set of all Gaussian states in (M) . We say that a unitary operator U 
inL^(R") is a Gaussian symmetry if, for any p £S„, the state C/p[/^ is also in Sn- All 
such Gaussian symmetries constitute a group ^^„. If a G C" and L S Sp{2n,'R) then 
the associated Weyl operator W(cic) and the unitary operator F(L) implementing the 
Bogolioubov automorphism of CCR corresponding to L are in (See Corollary 
3.5 in m.) The aim of this section is to show that any element U in is of the 
form XW {(X)r (L) where A is a complex scalar of modulus unity, (X E C" and L e 
Sp{2n,M.). This settles a question raised in |4j|. 

We begin with a result on a special Gaussian state. 

Theorem 6. Let si > S2 > ■ ■ ■ > s„ > be irrational numbers which are linearly 
independent over the field Q of rationals and let 

" + 

p. = p,(0,0,5) = n(l-e-v)e .^1 

7=1 

be the Gaussian state in L^(R") with zero position and momentum mean vectors and 
covariance matrix 

D = diag{di,d2,...,dn) 

with dj = jCOthjSj. Then a unitary operator U in L'(M") has the property that 
is a Gaussian state if and only if, for some OC G C", L G Sp{2n,R) and a 
complex-valued function ji of modulus unity on Z'J_ 

U = W(a)r(L))3(fl|fli,4fl2, . . . ,fllfl«) (4.1) 
where Z+ = {0, 1,2, . . .}. 

Proof. Sufficiency is immediate from Corollary 3.3 and Corollary 3.5 of H. To 
prove necessity assume that 



D 
D 



UpsU^ =Pg{e,m,S') 



(4.2) 
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Since a in L^(M) has spectrum Z+ and each eigenvalue k has multiplicity one ||2l 

" t 

it follows that the selfadjoint positive operator ^ SjUjaj, being a sum of commuting 



self adjoint operators Sja-aj, I < j <n has spectrum < ^ sjA:; 

;=i 



/t; € Z+ V ^ with 



each eigenvalue of multiplicity one thanks to the assumption on {sj, 1 < j < n}. 
Since Ps and UpsU^^ have the same set of eigenvalues and same multiplicities it 
follows from Theorem|4]that 



upsU-^ =w{z)r{M)-^p,r{M)w{zy^ 

where z € C", M e Sp{2n,R), t = (fi,f2, . . .,t„Y and 



(4.3) 



" + 

0, = Y{{\-e-'ne . 



Since the maximum eigenvalues of Ps and p, are same it follows that 



')• 



Since the spectra of pj and p, are same it follows that 



Choosing = (0,0, ... ,0, 1,0, . . . ,0)^ with 1 in the k-th position we conclude the 
existence of matrices A,B of order nxn and entries in Z+ such that 

r=As, s=Bt 

so that fiAs ~ s. The rationally linear independence of the i/s implies BA — I. This 
is possible only if A and B = A^' are both permutation matrices. 
Putting V = r{M)W{z)'^U we have from (|43]l 

Vps = p,V. 

Denote by |^) the vector satisfying 

QjOj \k) — kj \k) 
where |*;) = \ki)\k2) ■ ■ ■ \k„). Then 

Vps \k) = Y\{\-e-'y)e-^''''i V \k) 

.7=1 

= p,y|fe), fcez';. 
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Thus V\k) is an eigenvector for p, corresponding to the eigenvalue 

j=i 

= n(l-e-0)e-''^*. 

Hence there exists a scalar j3 {k) of modulus unity such that 

V\k) ^Pik) \Ak) 

= r(A)j3(fl|ai,4a2, • • ■,ala„) \\k) V A; G Z';. 

where F (A) is the second quantization of the permutation unitary matrix A acting in 
C". Thus 

U = W(z)r(M)^r(A)j3 (fl|fli,4fl2, . . .,ala,). 
which completes the proof. □ 

Theorem 7. A unitary operator U in l2(R«) is a Gaussian symmetry if and only if 
there exist a scalar A of modulus unity, a vector CC in C" and a symplectic matrix 
L g 5/?(2n,M) such that 

u = Xw{a)r{L) 

where W(OC) is the Weyl operator associated with CC andr{L) is a unitary operator 
implementing the Bogolioubov automorphism of CCR corresponding to L. 

Proof. The if part is already contained in Corollary 3.3 and Corollary 3.5 of Q. 
In order to prove the only if part we may, in view of Theorem|6l assume that U = 
P [a\a\,a\a2^ ■ ■ ■ ,a\an) where j3 is a function of modulus unity on U\. If such a U 
is a Gaussian symmetry then, for any pure Gaussian state | y/) , f/ 1 V^) is also a pure 
Gaussian state. We choose 

W)^e-^M'\e{u))=W{u)\e{0)) 

where m = (mi , M2, . . . , Wn)"^ G C" with uj 7^ V y. By our assumption 

1^'^) = e-ill"ll'j3(fl|fli,ajfl2, ■ ■ . ,ala„)\e{u)) (4.4) 

is also a pure Gaussian state. By Corollary|2] 3a G C", a unitary matrix A of order 
n and Ay > 0, 1 < y < n such that 

|v/'>=W(a)^WK)kAj---la„)- (4.5) 

Using ( 14.4b and ( I4.5l l we shall evaluate the function/(z) = (\//'|e(z)) in two different 
ways. From ( 14.4b we have 
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l|l.,l|2 , 
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/(z) =e (e(M) j3(fl|ai,4fl2, e(z)) 

= ^-^11""' L %i^^^(«-iZi)'^'...(«-„z„)'"l^i^2---U (4.6) 



Where |^i^2---^«> = |fci>|^2>---|^«> and |e(z))= I ^|^)forz€C. 

ytGZ+ 

Since |j3(jk)| = 1, ( |43] l implies 



|/(z)|<exp|-i||«||2+f 



(4.7) 



From the definition of \ex) in Corollary |2] and the exponential vector \e{z)) in 
L^fR) one has 
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I fx'- -I 



exp 



2 V A2 + 1 



A>0, 



This together with ( |4.5l l implies 

/(z) = (a,®a2®---®a„ |r(A-')w(-a)e(z)) 

= e<«l^>-5ll«ll'(e;Lj ® ^2 ® ■ ■ ■ ® a„ I e{A-\z+a))) 



which is a nonzero scalar multiple of the exponential of a polynomial of degree 2 
in zi ,Z2, • • • ,Zn except when all the A/s are equal to unity. This would contradict the 
inequality (14.61 1 except when A/ = 1 V j. Thus A/ = 1 V j and ( 14.5b reduces to 

I v/'> -w(a)r(A) |e(0) 
= .-^ll«in.(a)). 



Now (I4.4l l impUes 



J3(fl|fli,fl2a2, ■ ■ ■,alan) |e(M)) 



or 



I5{ki,k2,...,kn)\kik2...kn) 



/kil---Vk„l 



Thus 
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\U\ J \l'n 

Since |j3 (*) | = 1 and uj ^ OV j it follows that 1^1 = 1 and 

where 0/s are real. Thus j3(ajfli ,0202, ■ ■ ■ —^{D), the second quantiza- 

tion of the diagonal unitary matrix D = diag(e'^i , e'^ , . . . , e'^" ) . This completes the 
proof. □ 
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